How Long Could We Live? 

Ilja Dorsner^ and Pavel Fileviez Perez 1 ' 2 ' 2 ^ 

1 The Abdus Salam International Centre for Theoretical Physics 
Strada Costiera 11, 34014 Trieste, Italy. 
2 Pontificia Universidad Catolica de Chile. Facultad de Fisica, Casilla 306 

Santiago 22, Chile. 
3 Centro de Fisica Teorica de Particulas. Departamento de 
| Fisica. Instituto Superior Tecnico. Avenida Rovisco Pais,l 

1049-001 Lisboa, Portugal. 

(N 

We investigate model independent upper bounds on total proton lifetime in the con- 
text of Grand Unified Theories with the Standard Model matter content. We find them 

7— I 

(N ! to be r p < 1.5±g:l x W 39 (MX 'XJ CY)1 (0.003 GeV 3 /«) 2 years and r p < 7.lt°;° x 

"H- ' 10 36 (Mx/1 f 6 GcV)4 (0.003 GeV 3 /a) 2 years in the Majorana and Dirac neutrino case, re- 

a GUT 

qq | spectively. These bounds, in conjunction with experimental limits, put lower limit on the 

mass Mx of gauge bosons responsible for the proton and bound-neutron decay processes. 
For central values of relevant input parameters we obtain Mx > 4.3 x 10 14 y/acuT GeV. Our 
result implies that a large class of non-supersymmetric Grand Unified models, with typical 

P 1 

^5 values cxgut ~ 1/39, still satisfies experimental constraints on proton lifetime. Our result 

Oh! is independent on any CP violating phase and the only significant source of uncertainty is 

Oh, associated with imprecise knowledge of a — the nucleon decay matrix element. 

> 

jj ■ I. INTRODUCTION 

Grand Unified Theories Q, Q, Q, 0] (GUTs) are the most appealing extensions of the Standard 
Model (SM) of strong and electroweak interactions. Being founded on the ideas of force and matter 
unification they always generate two predictions regardless of their exact realization; one is the 
gauge coupling unification and the other is the proton decay. Of the two it is the latter that offers 
the only unambiguous way to test GUTs However, despite systematic experimental search it 
has not been observed so far [f], Q, Q]. Even if it is observed, a clear test of GUT might prove 
difficult due to inherent model dependence of all relevant proton decay contributions 
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Regardless of that, it is worth asking whether we can expect the test of the GUT idea through 
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proton decay experiments with certainty. 

There are several generic contributions to nucleon decay in GUTs. (For an incom plet e list of 
various studies on proton decay constraints on different unifying theories see 



20, 



21 



22 



23j].) In the non-super symmetric case the most important ones are the Higgs and gauge 



d = 6 contributions. In super symmetric theories there are two more contributions that generically 
predict too rapid proton decay. These are the d = 4 and d = 5 contributions. Of course, the four 
contributions we mention do not encompass all the possibilities. (For example, presence of extra 
Higgs representations such as 15 or 10 in an SU (5) GUT can lead to additional contributions 
through mixing of ap pro priate components of these representations with the triplet partners of the 



usual Higgs doublets 



241 . |25| | . It is also possible to have sizable contributions without any reference 



to the GUTs if the theory is super symmetric [26(.) But, they are certainly the most generic ones. 

It may come as a surprise that despite their multiplicity and diverse origins all of these con- 
tributions can in principle be completely suppressed or forbidden except the gauge d = 6 ones. 
For example, the so-called matter parity forbids the dangerous d = 4 contributions and there are 
numerous different ways to efficiently suppress the d = 5 operators and Higgs d = 6 operators in 
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3. 



realistic scenarios. (For discussion on suppression of d = 5 operators see for example 
In essence, the most promising tests of GUTs can be done through the gauge d = 6 contributions. 

The idea of using the gauge d = 6 dominated branching ratios for the two-body nucleon decays to 
distinguish between different GUT models of fermion mass has been introduced in the pioneering 
work of De Rujula, Georgi and Glashow j^]. Their idea has been revisited and elaborated on 
more recently. Namely, it has been shown that it is possible to make clear test of an y G UT with 
symmetric Yukawa couplings through the nucleon decay channels into antineutrinos |28j |. Similar 
conclusions j^] also hold in the context of flipped SU(5) j^ . 30 . 31, 32]. There, the clear test 
requires symmetric Yukawa couplings in the down-quark sector only. (Flipped SU(5) is to be 
considered a true GUT in the case of further embedding in S'O(IO).) 

But, in general, even the gauge d = 6 contributions can be significantly suppressed if not set to 



zero. For example, one can completely rotate them away in the flipped SU(5) context 33]. The 
relevant contributions there, which we refer to as the "flipped SU(5) contributions" for obvious 
reason, represent only one half of all possible gauge d = 6 contributions in GUTs. The other half, 
which we refer to as the ll SU (5) contributions" which is due to exchange of proton decay mediating 
gauge fields present in an SU(5) gauge group cannot be rotated away without a conflict with the 



measurements on fermion mixing 



341 ] . Nevertheless, it is worth investigating how efficiently one 



can suppress these contributions, too. Since there are no other gauge d = 6 contributions besides 
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the two we mention, this allows us to set an absolute upper bound on nucleon decay lifetimes. 
Crucial importance of those bounds lies in the fact that they are the only way to know if there is 
ever hope to test the idea of grand unification with certainty through proton decay experiments. 
Even if these bounds turn out to be beyond the experimental reach they set correct lower limit on 
Mqut through an absolute lower bound on the mass of the nucleon decay mediating gauge bosons. 
In other words, they are the bounds that can tell us which GUT scenarios are a priori ruled out 
by experimental data. In what follows we concentrate on GUTs with the SM matter content, i.e., 
the three generation case, due to their phenomenological relevance. 



II. LOOKING FOR AN UPPER BOUND ON THE TOTAL PROTON LIFETIME 

To establish an upper bound on the total proton lifetime we first critically analyze all possi- 
ble gauge d = 6 operators contributing to proton decay. Again, we concentrate solely on these 
contributions since all other contributions can be set to zero. 

Proton lifetime induced by superheavy gauge boson exchange can be written as follows 

r p = C M x a ~ 2 UT m- 5 , (1) 

where C is a coefficient which contains all information about the flavor structure of the theory. 
Mx is the mass of the superheavy gauge bosons, ugut = 9gutI^^ wnere 9GUT is the coupling 
defined at the GUT scale (the scale of gauge unification). To find a true upper bound on the total 
lifetime we need to find the maximal value for the C coefficient. Then, for a given value of Mx 
and olgut we can bound the GUT scenario prediction for the nucleon lifetime. 

The relevant gauge d = 6 operators contributing to the decay of the proton, in the physical 
basis 28], are: 

°{ea,dp) = c(e%,dp) e ijk uf ^ uj e% 7 M d k p, (2a) 
0(e a ,dp) = c(e a ,dp) e ijk uf 7 M uj d^ 7^ e a , (2b) 
0(vi,d a ,dp) = c(vi,d a ,dp) e ijk uf 7^ d ja d^ 7^ v h (2c) 
0(y x ,d a ,dp) = c(vf,d a ,dp) e ijk dfp 7^ Uj vf 7^ d ka , (2d) 
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where the relevant coefficients are given by: 



c(e a , dp) 
c(e a ,dp) 



la 



(3b) 



(3a) 



cfa,d a ,c$) = f^{ViV UD ) la {V z V EN f l + kiv£ a {ViV UD V}V 3 V EN ) 11 , a = 1 or (3 = 1, (3c) 
c(vf,d a ,dC) = kl[{V^ D f\U EN V 2 ) la + vt(U j EN V 2 vi D )% a = 1 or (3 = 1. (3d) 

The mixing matrices Vi = utu, V 2 = E ] C D, V 3 = DlE, V 4 = Dt,D, V UD = U^D, V EN = 



E^N, and Uen = EqNq- a, (3 = 1,2, I = 1,2,3, while and k are the color indices. (Our 
convention for the diagonalization of the up, down and charged lepton Yukawa matrices is specified 



by XfgYuU = Y^ e , D T C Y D D = Y^ g , and E%Y E E = Y^ iag .) The quark mixing is given by V UD = 



WD = K\VckmK2, where K\ and K 2 are diagonal matrices containing three and two phases, 
respectively. The leptonic mixing Ven = K^VfK^ in case of Dirac neutrino, or V E n = K^V^ 1 in 
the Majorana case. Vf and V^ M are the leptonic mixing matrices at low scale in the Dirac and 
Majorana case, respectively. The gauge d = 6 operators have to be run from the GUT scale down 
to 1 GeV, i.e., the proton decay scale, and the appropriate amplitude computed in the usual way. 



(For details, see for example |35|.) 

In the above expressions k\ = gGUTM7 x Y y and k 2 = gci/rM^i Y ')i where M^ X) y), M(X',Y') 
~ Mqut are the masses of the superheavy gauge bosons. All terms proportional to k\ are obtained 
when we integrate out (X, Y) = (3, 2, 5/3), where X and Y fields have electric charge 4/3 and 1/3, 
respectively. These are the fields appearing in theories based on the 577(5) gauge group. Thus, 
we call their contributions the "577(5) contributions". Integrating out (X',Y') = (3,2,-1/3) 
we obtain the terms proportional to k 2 . These contributions we refer to as the "flipped 577(5) 
contributions" since they appear in the flipped SU(5) scenario. The electric charge of Y' is —2/3, 
while X' has the same charge as Y. Again, there are no other gauge contributions in any GUT 
besides these. 

Minimization of the total decay rate represents formidable task since there are in principle 42 
unknown parameters. To face the challenge we look for a solution where the "5C/(5) contributions" 
and "flipped SU(5) contributions" are suppressed (minimized) independently. Since we expect that 
in general the associated gauge bosons and couplings have different values this is also the most 
natural way to look for the minimal decay rate value. Moreover, the bounds obtained is such a 
manner will be independent of the underlying gauge symmetry. 
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The "nipped SU(5) contributions" are set to zero by the following two conditions |33{ |: 

Vf a = {D ] c Df a = 0, a = 1 or f3 = 1, (Condition I) 
{U ] c E) la = 0. (Condition II) 

(Condition I cannot be satisfied in the case of symmetric down quark Yukawa couplings.) Therefore, 
in the presence of all gauge d = 6 contributions, in the Major ana neutrino case, there only remain 
the contributions appearing in SU(5) models. These, however, cannot be set to zero j^] in the case 
of three generations of matter fields. But, as we now show, they can be significantly suppressed. 
There are two major scenarios to be considered that defer by the way proton decays: 

• There are no decays into the meson-charged antilepton pairs 

All contributions to the decay of the proton into charged antileptons and a meson can be 
set to zero . Namely, after we implement Conditions I and II, we can set to zero Eq. (JHEl by 
choosing 

Vi 1 = {U ] C U) 11 = 0. (Condition III) (4) 

(This condition cannot be implemented in the case of symmetric up-quark Yukawa cou- 
plings.) On the other hand, Eq. (|3a|) can be set to zero only if we impose 

(V 2 V^ D ) al = {E ] c U) al = 0. (Condition IV) (5) 

Therefore with Conditions I-IV there are only decays into antineutrinos and, in the Majorana 
neutrino case, the only non-zero coefficients are: 

c(u h d a ,d§) = kj (V^uo^iVsVENf 1 . (6) 

So, indeed, there exists a large class of models for fermion masses where there are no decays 
into a meson and charged antileptons. 

Up to this point all conditions we impose are consistent with the unitarity constraint and 
experimental data on fermion mixing. (In the SU(5) case we have to impose Conditions III 
and IV only.) We now proceed and investigate the decay channels with antineutrinos. From 
Eq. © we see that it is not possible to set to zero all decays since the factor {V\VuD) la can 
be set to zero for only one value of a in order to satisfy the unitarity constraint. Therefore 
we have to compare the following two cases: 
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Case a) (VxVud) 11 = (Condition V). 

In this case the chiral langragian technique yields: 



r a (p- 




= o, 


r a ( P - 


K + 9) 


= C(p,K) 


r a (n - 


* vr°F.) 


= o, 


r a (n- 


*K°v) 


= C(n,K) 


r a (n 


-> Wi) 


= o, 



1 + 



m r 



3m.B 



in 



i 2 



(D + 3F) 



l^CKM^CK 
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'CKM v CKM V CKM V CKM\ 



,2 i IT/21 12 
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3ms 



2 It/32 T/21 _ T/31 T/ 



22 1 2 
CKTM V CKM\ 



IT/31 12, IT/21 12 
I V CKM\ ^ I "OKTAfl 



where: 



C(a,b)= {m }_jjf A\ \a\ 2 fc?. 



(7) 
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Case b) (yiVc/D) 12 = (Condition VI). 

All the decays channels into antineutrinos are non-zero in this case. Associated decay 
rates are: 

IT/32 t/21 



T b (p -> 7T + Z/) 
T fe (p -> JT+i/) 
T 6 (n -> 7T°u) 

T b (n -> 
T fe (n -> 7/17) 



C(p,7r) [i + jD + f] 2 ^ m ^ m 



T/31 t/22 I 2 



2m p 
3m b 



IT/22 12 _i_ |T/32 12 
2 it/32 t/21 _ t/31 T/22 12 
IT/22 |2 i |T/32 |2 : 



C( n , tt) [1 + jP +F]2 r(p -> 7T+zy), 



C(n,iT) 



1 + — — (D-3F) 



2 it/32 t/ 21 _ T/31 



v 



22 



CK M v CKM v CKM v CKM 
IT/22 |2 i IT/32 |2 
\ V CKM\ ~ r \ v CKM\ 



C(n,rj) 



[1 + D-3F] 2 \V C 



2 It/32 t/21 t/31 T /22 1 2 



CKM V CKM " CKM* CKM\ 
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IT/22 12 I |T/32 |2 

I ckm\ \ v ckm\ 



The nice thing about these results is that they are completely independent of all CP violating 
phases including those of Vckm and Vi and any mixing angles beyond the CKM ones. (This 
is completely unexpected since there are in principle 42 different angles and phases that could 
a priori enter our results.) Also, in the limit VMsr ~~ all decay rates vanish as required in 



the case of three generations of matter fields 34]. To demonstrate these two properties we 
adopt the so-called "standard" parametrization of Vckm that utilizes angles 

^12, #23, #13, and a phase £13. (For example, in that parametrization Vq 3 km = e~ lSvi S13.) The 
relevant terms read V 32 KM V 21 KM -V C 1 KM V 22 KM = e iS ^s 13 , \V% 2 KM \ 2 + \V$ 2 KM \ 2 = c 2 12 + s 2 12 s 2 13 
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and I^c^jv/I 2 + Wckm\ 2 = s i2 + c i2 s i3> where = cos 6. L j and Sjj = sin#jj. Hence, all one 
needs to know are angles # 12 and #13 . 

We present numerical values of all relevant two body decay lifetimes for proton and bounded 
neutron decays in Tables HI and ITU respectively. Clearly, it is Case b) that gives the lowest 
total decay rate in the Majorana neutrino case. (We also include the Dirac neutrino case 
for completeness.) Lifetimes are given in units of / a GUT> wnere the gauge boson mass 
is taken to be 10 16 GeV. To generate these values we use m p = 938.3 MeV, D = 0.81, 
F = 0.44, m B = 1150 MeV, f w = 139 MeV, A L = 1.43, and the most conservative value 
a = 0.003 GeV 3 j^jj. Indicated uncertainties reflect the errors in measurement of angles 
9\2 and #13 only. These are well-known and their sines are: S12 = 0.2243 ± 0.0016, and 
S13 = 0.0037 db 0.0005 |4lj |. Note that the most poorly known parameter is actually a; the 
most recent QCD lattice calculations j^ . indicate that its value could be three times 
bigger than the value we use. If that result persists it would reduce the lifetime bounds we 
present by a factor of ten. 

TABLE I: Proton lifetimes in years for Majorana and Dirac neutrinos in units of Mj^ / 'a.Q UT , where the 
mass of gauge bosons is taken to be 10 16 GeV. 







Majorana 




Dirac 


Channel 


Case a) 


Case b) 


Case a) 


Case b) 


p — > 7T + P 


00 


ill x 1038 


5.4^1:1 x 1q38 


2.6t° Q 9 e x 10 38 


p — ► K + 9 


1.0±g;| x 10 38 


2.5t° Q 9 6 x 10 40 


6.8±°;g x 10 36 


7.2±°;° x 10 36 


TOTAL 


l.0t°j x 10 38 


5 *°-i:I x 1038 


6.7±H x 10 36 


7.lt°;° x 10 36 



TABLE II: Lifetimes for bounded neutrons in years for Majorana and Dirac neutrinos in units of M^/oq 
where the mass of gauge bosons is taken to be 10 16 GeV. 







Majorana 




Dirac 


Channel 


Case a) 


Case b) 


Case a) 


Case b) 


n — > ir°i> 


00 


1.0±°; 3 x 10 39 


l.l^; 4 x 10 39 


5.2t\i x 10 38 


n-> K°v 


1.1±°; 4 x 10 38 


6.7±i J x 10 39 


1.9±g;g x 10 36 


1.9l 9 ;° x 10 36 


n ^ r\v 


00 


1.5±g;| x 10 41 


1.6t 9 ; 5 3 x 10 41 


7.6l?;? x 10 40 


TOTAL 


l.lt°; 4 x 10 38 


8.8t 2 2 ; 9 x 10 38 


1.9t 9 ;° x 10 36 


1.9t° ° x 10 36 
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• There are no decays into the meson-antineutrino pair in the Majorana neutrino 
case 

Let us show that it is also possible to set to zero all nucleon decay channels into a meson 
and antineutrinos. After Conditions I and II, we could impose (ViVud) = (Condition 

VII) instead of Uj 11 = 0. (Again, these two equalities are exclusive in the case V^ M ^ 
0.) Therefore, in the Majorana neutrino case, there are no decays into antineutrinos (see 
Eq. (JHcjl ). In this case the property that the gauge contributions vanish as |V r ( i'|- M | — > is 
obvious since \V^\ = \ V^ KM \. We have to further investigate all possible values of and 
Vg a . Now, we can choose = and V3 a = 0, except for the case a = (3 = 2 (Condition 

VIII) . In that case there are only decays into a strange mesons and muons. Let us call this 
Case c). To understand which case gives us an upper bound on the total proton decay 
lifetime in the Majorana neutrino case, we compare the predictions coming from the Case 
b) and Case c). The ratio between the relevant decay rates is given by: 

T b (p^n+u) ~ 2(Cl2 + Sl2Sl3) (mj-miy [1 + D + F] 2 ~ ^ (8) 
Thus, the upper bound on the proton lifetime in the case of Majorana neutrinos indeed 
corresponds to the total lifetime of Case c). We find it to be: 

r P < 1.5±8;| x 10 39 (Mx/1 f 6QeV)4 (0.003 GeV 3 /a) 2 years, (9) 

a GUT 

where the gauge boson mass is given in units of 10 16 GeV. We explicitly indicate the depen- 
dence of our results on the nucleon decay matrix element. These bounds are applicable to 
any GUT regardless whether the scenario is supersymmetric or not. If the theory is based 
on SU(5) the above bounds are obtained by imposing Conditions VII and VIII. If the theory 
contains both SU(5) and flipped SU(5) contributions, in addition to these, one needs to 
impose Conditions I and II. 

We plot the proton bounds in the Mx-Q-gut plane for the Majorana (Dirac) neutrino case in 
Fig. ^ (j2J). Our results, in conjunction with the experimental limits on nucleon lifetime, set an 
absolute lower bound on mass of superheavy gauge bosons. Since their mass is identified with the 
unification scale after the threshold corrections are incorporated in the running 44| this also sets 



the lower bound on the unification scale. Using the most stringent limit on partial proton lifetime 



{jp > 50 x 10 32 years) for the p — ► ir°e + channel 41] and setting a = 0.003 GeV 3 [40j, we obtain 



M x > 4.3l°; 3 x 10 14 V^c^GeV, (10) 
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where cxgut usually varies from 1/40 for non-supersymmetric theories to 1/24 for supersymmetric 
theories. For example, if we take a non-supersymmetric value olgut = 1/39, we obtain 

M x >7x 10 13 GeV. (11) 

Again, this result implies that any non-supersymmetric theory with acuT = 1/39 is eliminated 
if its unifying scale is bellow 7.0 x 10 13 GeV regardless of the exact form of the Yukawa sector of 
the theory. Note that majority of non-supersymmetric extensions of the Georgi-Glashow SU(5) 
model yield GUT scale which is slightly above 10 GeV. Hence, as far as the experimental limits 
on proton decay are concerned, these extensions still represent viable scenarios of models beyond 
the SM. Region of Mx excluded by the experimental result is also shown in Figs. ^ and [21 
At this point the following two observations are in order: 

1. All three cases (Case a)-c)) yield comparable lifetimes (within a factor of ten) even though 
they significantly defer in decay pattern predictions; 

2. We use the most stringent experimental limit on partial proton lifetime as if it represents 
the limit on the total proton lifetime. Even though this is not correct (see discussion in j^l) 
it certainly yields the most conservative bound on Mx- 



One can easily extend our results to a class of orbifold GUT theories 



45 



where all matter 



fields live on an "unbroken" brane. In essence, to obtain the lower limit on the gauge boson 



mass in those theories, it suffices to multiply the limit presented in Eq. IjllJI by wtt/2. (This 
factor accounts for the fact that the two-body decay of the proton is due to exchange of an entire 
Kaluza-Klein (KK) tower of states j^] associated with the proton decay mediating gauge boson.) 
The bound obtained in such a way then corresponds to the limit on the compactification scale of 
extra dimension(s). (Recall that in orbifold GUTs the gauge bosons responsible for proton decay 
belong to the KK tower where the lightest gauge boson in the tower has the mass equal to the 
orbifold compactification scale.) Curiously enough, exact unification of gauge couplings in the 
five-dimensional S 1 /(Z2 x Zo)-type orbifold models usually requires the compactification scale to 



be slightly above 10 J ! GeV" 18. |9j,Q. This would imply that the orbifold GUT theories with the 
matter fields all located on the "unbroken" brane could soon be completely ruled out if the proton 
decay is not observed in the next generation of the proton decay experiments. 

In order to complete our analysis let us finally demonstrate the possibility to set to zero the 
Higgs d = 6 and d = 5 contributions. The triplets T = (3, 1, —2/3) and T = (3, 1,2/3) have the 
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FIG. 1: Isoplot for the upper bounds on the total proton lifetime in years in the Majorana neutrino case 
in the Mx-olgut plane. The value of the unifying coupling constant is varied from 1/60 to 1/10. The 
conventional values for Mx and olgut in SUSY GUTs are marked in thick lines. Experimentally excluded 
region is given in black. 

following interactions: 

W T = j d 2 e[[QAQ + U c CE C + D c E N c ] f + [ Q B L + U c D D c } f } + h.c. 

(12) 

Choosing = —A^ and = 0, except for i = j = 3, the Higgs d = 6 and d = 5 contributions 
are indeed set to zero. It is also possible to have SUSY scenarios where the d=5 operators are 
strongly suppressed by particular realization of superparticle spectrum In any case, even if 
SUSY is realized at low energies we are sure that the upper bound is coming from the gauge d = 6 
contributions. 
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FIG. 2: Isoplot for the upper bounds on the total proton lifetime in years in the Dirac neutrino case in the 
Mx-acuT plane. The value of the unifying coupling constant is varied from 1/60 to 1/10. The conventional 
values for Mx and aai/T in SUSY GUTs are marked in thick lines. Experimentally excluded region is given 
in black. 

III. SUMMARY 

We have investigated the possibility of finding an upper bound on the total nucleon de- 
cay lifetime in the context of grand unified theories with the Standard Model matter content. 
This bound originates from the gauge d = 6 contributions, since all other contributions are 
quite model dependent and can always be suppressed. In the Majorana neutrino case the 
bound is t p < 1.5 x 1 39 (^Af GeV ) 4 (0,003 GeV 3 /«) 

years, while in the Dirac neutrino case 

16 4 aGUT 

T p < 7.1 x 10 36 (M*/ 1 ® — GeV) QQ3 QeV 3 /a) 2 years. These bounds are valid in both supersymmet- 

a GUT 

ric and non-supersymmetric scenarios and are grand unifying gauge group independent. Moreover, 
there is no dependence of our results on CP violating phases nor any angles beyond those of CKM. 
Our bounds are very useful for two reasons. Firstly, in the context of realistic grand unified theories 
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they indicate whether it is possible to test these theories in their entire flavor parameter space with 
certainty through proton decay experiments. Secondly, they put an absolute lower bound on the 
mass of proton decay mediating gauge bosons. We obtain Mx > 4.3^ J x 10 14 ^acuT GeV for a 
reasonable set of input parameters. Since this mass is usually identified with the unifying scale 
through threshold matching conditions our bounds can be interpreted as the lower bounds on the 
GUT scale itself. We have also addressed implications our bounds have on the popular class of the 
so-called "orbifold" models. 
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